Quark-hadron phase transition in a neutron star under strong magnetic fields 
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We study the effect of a strong magnetic field on the properties of neutron stars with a quark- 
hadron phase transition. It is shown that the magnetic field prevents the appearance of a quark 
phase, enhances the leptonic fraction, decreases the baryonic density extension of the mixed phase 
and stiffens the total equation of state, including both the stellar matter and the magnetic field 
^\ ' contributions. Two parametrisations of a density dependent static magnetic field, increasing, re- 

, spectively, fast and slowly with the density and reaching 2 — 4 x 10^* G in the center of the star, 

' are considered. The compact stars with strong magnetic fields have maximum mass configurations 

' with larger masses and radius and smaller quark fractions. The parametrisation of the magnetic 

field with density has a strong infiuence on the star properties. 
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PACS numbers: 26.60.+C, 12.39.Ba, 21.65.+f, 97.60. Jd 



I. INTRODUCTION 

y ; Neutron stars with very strong magnetic fields of the order of 10"'^'' — 10^^ G are known as magnetars and they are 
■ ■ believed to be the sources of the intense gamma and X rays [1-3]. If the magnetic field at the center of the star has 
the same magnitude as the one measured at the surface of a magnetar, then the influence of the magnetic field on 
the star's mass and radius is negligible. However, according to the scalar virial theorem the maximum interior field 
strength could be as large as 10^* G for a compact star [4]. For a field of this magnitude or larger considerable effects 
will be induced on EOS. 

It has been shown by several authors that the magnetic fields larger than B = 5 x 10^* G will affect the EOS of 
compact stars [5, 6]. In particular field-theoretical descriptions based on the non-linear Walecka model (NLWM) [7] 
were used. The generalization of these works to density-dependent hadronic models [8] or the inclusion of the scalar- 
isovector 5 mesons [9] have also been applied to the description of hadronic stars with intense magnetic fields. Quark 
^) ■ stars with strong magnetic fields have been described within the MIT bag model [5, 10, 11] and the SU(2) version 
^\ of the Nambu-Jona-Lasinio model (NJL) [12], an effective model which includes the chiral symmetry, a symmetry of 
O QCD. 

l-^ Compact stars are complex systems which may contain in their core a pure quark phase or a non-homogenous 

mixed quark-hadron phase [13]. We would like to investigate whether the presence of a magnetic field could affect the 
presence of a quark core. The influence of strong magnetic fields on quark-hadron phase transition was first discussed 
in Ref. [14], using a Dirac-Hartree-Fock approach within a mean-field approximation to describe both the hadronic 
and quark phases. For the hadronic matter they took a system of protons, neutrons and electrons and for the quark 
phase the MIT bag model with one-gluon exchange. A very hard quark EOS was used so that the hybrid star did 
not have a quark core. They have concluded that compact stars have a smaller maximum mass in the presence of 
strong magnetic fields, a result that does not agree with other works where hadronic stars [15] or quark stars [12] in 
the presence of strong magnetic fields, have been studied. 

In the present work we will study the effect of a strong magnetic field on the deconfinement phase transition inside 
a compact star. We will consider two EOS for the quark phase, a softer and a harder one, chosen in such a way 
that the star will have at least a small quark core if no magnetic field exists. For the hadronic phase we consider 
a standard hadronic EOS, the GMl parametrisation of the non-linear Walecka model proposed in [16], and include 
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the baryonic octet. The largest magnetic field measured at the surface of a star is ^ 10^^ G but nothing is known 
about the magnitude of magnetic fields in the interior of a star. We will consider a magnetic field whose magnitude 
increases from a value of ~ 10^^ at the surface to a maximum value of ~ 4 x 10^^ G at the center. In particular, we 
want to analyse how the parametrisation of the field affects the star properties, such as the radius and mass. 

This work is organized as follows: we make a brief review of the formalism used for the hadronic, mixed and quark 
phases. Next we discuss the parametrisation used for the magnetic field and present and discuss the results obtained 
for the equation of state of stellar matter and the mass/radius properties of the corresponding compact star families, 
for several values of the maximum magnetic field. At the end we will draw some conclusions. 



II. HADRON MATTER EQUATIONS OF STATE 

For the description of the EOS of neutron star matter, we employ a field-theoretical approach in which the baryons 
interact via the exchange of a — uj — p mesons in the presence of a static magnetic field B along the 2;-axis. The 
Lagrangian density of the NLW model (GMl) can be written as 

jC = J2^b + jCm + Y.Cl. (1) 

b I 

The baryons, leptons {l=e, fi), and mesons (cr, uj and p) Lagrangians are given by 

-Cb = {iliid^ - qtliiA'^ -mb + QabO^ - Qwbln^^ - 9pbT3tlf^P'' - ^P-Nl^b^^l^iyP""^ *6 

- ^F^-^F„. + ^mlp,p^-^P>-P,, (2) 

where ^E*;, and ipi are the baryon and lepton Dirac fields, respectively. The index b runs over the eight lightest baryons 
n, p, A, S~, T,^, S+, E~ and E^. (Neglecting the baryonic decuplet including the fl~ and the A quartet, which appear 
only at very high densities, does not affect our conclusions). The baryon mass and isospin projection arc denoted 
by TUb and t^^, respectively. The mcsonic and electromagnetic field tensors are given by their usual expressions: 
^tiv = dfiOJi, — d,yUJf_i, Pfii, = dfj,pi, — di,p^, and F^i, = d^A^ — dvA^. The baryon anomalous magnetic moments (AMM) 
are introduced via the coupling of the baryons to the electromagnetic field tensor with a^v = f [7/i) lu] and strength 
Kb- The electromagnetic field is assumed to be externally generated (and thus has no associated field equation), and 
only frozen-field configurations will be considered. The interaction couplings are denoted by g, the electromagnetic 
couplings by q and the baryons, mesons and leptons masses by m. The scalar self-interaction is taken to be of the 
form 

U {cr) = ^bnin {gaNcr f + {gaNcr f ■ (3) 

The parameters of the model are the nucleon mass M = 939 MeV, the masses of mesons m,^, m^^, nip and the coupling 
parameters. 

The mcson-hypcron couplings arc assumed to be fixed fractions of the meson- nucleon couplings, giH = XingiN, 
where for each meson i, the values of XiH are assumed equal for all hyperons H. The values of xm are chosen to 
reproduce the binding energy of the A at nuclear saturation as suggested in [16], and given in Table II. 

From the Lagrangian density in Eq. (1), we obtain the following meson field equations in the mean-field approxi- 
mation 

2 dU (a) „ „ 

m^cr + = 2^ 9abPb = 9aN 2^ X^bPb (4) 

6 6 

b b 

"^pP^ = ^9pbT3,Pb = 9pN^XpbT3^pl (6) 
5 6 

where cr = (cr) , = (a;°) and p = (p") are the nonvanishing expectation values of the mesons fields in uniform 
matter. 

The Dirac equations for baryons and leptons are, respectively, given by 
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TABLE I: Static properties of baryons considered in this study. The mass, charge of baryon and strange charge of species b 
are denoted by mt, Qb and g^, respectively. The baryonic magnetic moment is denoted by /^(,. The baryonic AMM is given by 



Kb = {nb/ UN — qbmp/mb), where fj,N is the nuclear magneton. 



b mb(McV) qb{e) g| Hb/HN Kb 

p 938.27 1 2.97 1.79 

n 939.56 -1.91 -1.91 

A'' 1115.7 -1 -0.61 -0.61 

S+ 1189.4 1 -1 2.46 1.67 

S° 1192.6 -1 1.61 1.61 

E" 1197.4 -1 -1 -1.16 -038 

E° 1314.8 -2 -1.25 -1.25 

E- 1321.3 -1 -2 -0.65 0.06 



[iT/.S" - 967/*^^ -ml -Jo {gu:<^° + QpT^y) - ^fiNKb'J^.uF'"'] = (7) 

ih^.d^-qn^,A^'-mi)i^i = (8) 

where the effective baryon masses are given by 

ml = mb- ga(j (9) 

and pI and pi are the scalar number density and the vector number density, respectively. For charge-neutral, neutrino 
free, /3-equilibrated matter, the following conditions are satisfied: 

Ma = Ms" = Mho = /x„. 
Me- = Mh- = l^n+f^e, 
IJ-IX = IJ-e, (10) 

Y.^bPl + Y.^lp'i=Q- (11) 

b I 

The energy spectra for charged baryons, neutral baryons and leptons (electrons and muons) are, respectively, given 

by 



El,s = y fcf + (^\fm^+Mq^ - spNUbB^ + g^bio° + Tz,9pbP° (12) 

El = '^kl+ (^K^ + k:i + kl- suNKbB^ + g^boj'^ + T^.gpbP^ (13) 
Ei^, = ^kl + m'i + 2v\qi\B (14) 

where v = n+\ — sgn{q)^ = 0, 1, 2, . . . enumerates the Landau levels (LL) of the fermions with electric charge q, the 
quantum number s is -1-1 for spin up and —1 for spin down cases. 

For the charged baryons, we introduce the effective mass under the effect of a magnetic field 

fhl = ^Jml^ + 2v\qb\B - sfj,NKbB (15) 
the expressions of the scalar and vector densities are, respectively, given by [6] 



^F,i/,s + Ep 
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27r2 



(16) 



i/=0 s 



where ^ is the Fermi momenta of charged baryon 6 with quantum numbers v and s. The Fermi energy is 
related to the Fermi momenta ^ by 



(fcks) =(i?^) -(mgf 



(17) 



For the neutral baryons, there is no LL quantum number so the Fermi momenta is denoted by fc^^, and the Fermi 
energy E\ is given by 



hi 



with 



mh = ml — sunhB. 
The scalar and vector densities of the neutral baryon b are, respectively, given by 



(18) 
(19) 



Pb = 

Pi = 



1 
2^ 



E 



i;^4,-m^ln 



4,. + E'p 



rub 



1 



{k^-^af ~ \siiiqHbB i nikkpa + (Ep)^ i arcsin 



3 V r,.j 2 

The vector density for leptons is given by 



27r2 



where kp ^ ^ is the lepton Fermi momenta, which is related to the Fermi energy Ep by 

{kF,u,s)^ = (Epf -fhf, l = e,ii, 



(20) 



(21) 



(22) 



with fhf = mf + 2i'\qi\B. The summation in v in the above expressions terminates at Vmax, the largest value of u for 
which the square of Fermi momenta of the particle is still positive and which corresponds to the closest integer from 
below defined by the ratio 



(4 



2\qi\B 



leptons 



{E^ + s^iN KbBf-mf 



2\qb\B 

The chemical potentials of baryons and leptons are defined as 



, charged baryons. 



Mb 



Ep + g^6t^° + gpbT3^p° 



E'f = J {k^F,.,s) +m^ + 2t^\qi\B. 



(23) 
(24) 



We solve the coupled Eqs. (4)-(ll) self-consistently at a given baryon density p = J2b Pb ™ presence of strong 
magnetic fields. The energy density of neutron star matter is given by 



(25) 



l=e,iJ, 
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where the energy densities of charged baryons eg, neutral baryons s^, and leptons e; have, respectively, the following 
forms 



sr. = 



\qb\B 
47r2 
1 

4^2 



EE 



1^=0 s 
1 „ 1 



Ic" _|_ 



i E fl^^. (^f)' - Is^^N^^B {E^,f ( arcsin ( § ) - | 



-sfiNKbB + -TTib ) I mbkpgEp + In 



Si 



\H\B 
4772 



EE 



kF.,sEF+ffitln 



i/=0 s 

The pressure of neutron star matter is given by 

= E ^'^^ - = E 



(26) 



(27) 



where the charge neutrality and /^-equilibrium conditions are used to get the last equality. The total energy density 
and the total pressure of the system are given, by adding the corresponding contribution of the magnetic field. 



^2 

s2 

P = Pr. + ^. 



(28) 



III. QUARK MATTER EQUATIONS OF STATE 

The MIT Bag model has been extensively used to describe quark matter. In its simplest form, the quarks are 
considered to be free inside a Bag and the thermodynamic properties are derived from the Fermi gas model. In the 
presence of a strong magnetic field Eqs. (28) are still valid with the energy density and pressure of quark stellar 
matter, and quark density, respectively, given by 



i=u,d,s,e 




z=u.d,s.e 



In 



In 



Hi + y 












Mi + ] 






(m«)' 



E5VMf-(M«)' 



+ Bag, 



Bag, 



(29) 



where M^*' = + 2i'\qi\B and v runs over the allowed LL, Qi denotes the degeneracy of the i-th species, i.e., six 
for the quarks and two for the leptons, the quarks masses and Bag represents the bag pressure. 

We use niu = ma = 5.5MeV, m,, = 150.0McV and two values for the Bag (165MeV)^ and (180MeV)^. In what 
follows, the baryon density p is defined as p = J2q Pql'^- 

In a star with quark matter we must impose both beta equilibrium and charge neutrality. For /3-equilibrium matter 
we must add the contribution of the leptons as free Fermi gases (electrons and muons) to the energy density, pressure 
and entropy density. The relations between the chemical potentials of the different particles are given by 



Ms = = + Me, Me = M/i- 



(30) 
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In terms of neutron and electric charge chemical potentials /i„ and /Xei one has 



12, 11 

/"« = ^IJ-n - ^Me, and IJ,g = IJ,d = -fj,n + -/ie- (31) 



For the charge neutrality we impose 



Pe + P/n = 2 (2p„ - pd- Ps) ■ (32) 



IV. MIXED PHASE AND HYBRID STAR PROPERTIES 



We study the hadron-quark phase transition which may occur in the core of massive neutron stars. We will use the 
Gibbs criteria and global charge conservation to describe the mixed phase of hadronic and quark matter [13], which 
for two conserved charges, the electric charge and the baryonic charge, are given by 

MHP,i = PQPa = Pi^ i = n,e, 

PapiPHp) = Pqp{pqp)- (33) 

where the subscripts HP and QP denote, respectively, the hadronic and the quark phase. Imposing the condition 
of global conservation of the electric and baryonic charges, i.e, both hadronic and quark matter are allowed to be 
separately charged and have different baryonic densities. We impose 

Xp'^^ + {1 - x) pf^ + p'c = 0, 

^pQP + (l_^)pHP ^ p (34) 

where x if' the volume fraction occupied by the quark phase and (1 — x) the volume fraction occupied by the hadron 
phase, pI (p*) arc the electric (baryonic) charge densities of quark and hadron phases and p^. is the lepton electric 
charge density. The total energy density in the mixed phase reads 

s = X£^^ + (l-x)£"^ + £;- (35) 

The mixed phase is charactrizcd by a value of x which varies between 0, at the onset of the mixed phase, and 1, at the 
onset of the quark matter. The baryonic densities pi and p2 given in Tables III and IV, correspond, respectively, to 
X = and x = 1 ^^nd define the onset of the mixed phase (pi) and the onset of the quark phase (^2)- The mixed phase 
is determined by simultaneously imposing the Gibbs conditions (33), charge conservation conditions (34), chemical 
equilibrium conditions (10), (11), (31) and (32), and the field equations (4-8). 

However, we should point out that the Gibbs construction is an approximation that although taking correctly into 
account the existence of two charge conserving conditions does not take into account surface effects and the Coulomb 
field [17, 18]. A complete treatment of the mixed phase requires the knowledge of the surface tension between the two 
phases which is not well established and may have a value between 10-100 McV/fm^ [18]. The Gibbs construction 
gives results close to the ones obtained with the lower value of the above surface tension range, and is recovered 
for a zero surface tension. We could have also considered the Maxwell construction, for which only baryon number 
conservation is considered, i. e. the first condition of Eq. (33) is only satisfied for i = n. It has been shown in [18] that 
the Maxwell description of the mixed phase gives a good description if the surface tension is very large. Therefore, 
for comparison we will also determine the mixed phase using a Maxwell construction, see Tables V and VI. 

The EOS for the mixed phase are then constructed. Consequently, we can compute the properties of the neutron 
star. In the present study we consider very strong magnetic fields which have a contribution to the total pressure at 
least as large as the contribution coming from matter. This means that stars with non-spherical configurations should 
be considered [15, 19]. We will, however, consider spherical stars and use the Tolman-Oppenheimer-Volkoff (TOV) 
equations for the structure of a static, spherically symmetric star to obtain the mass-radius (M-R) relation of neutron 
star. With this approximation our results and, therefore, conclusions should be taken with care. 



V. RESULTS AND DISCUSSION 



In order to study the effect of strong magnetic fields on the structure of neutron star we use the GMl parametrisation 

of the NLW model given in Table. II, [16] . Since, to date, there is no information available on the interior magnetic field 
of the star, we will assume that the magnetic field is baryon density-dependent as suggested by Ref. [14]. The variation 
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TABLE II: The parameter set GMl [16]used in the calculation. 



Po 


-B/A 




gaN/n-lc 




gpN/mp 














(MeV) 


M*/M 


(fm) 


(fm) 


(fm) 






XpH 


b 


c 


0.153 


16.30 


0.70 


3.434 


2.674 


2.100 


0.600 


0.653 


0.600 


0.002947 


-0.001070 




FIG. 1: (Color online). The magnetic field B defined in (36) as function of the baryonic density for several values of the 
parameter Bo and for the slow (/3 = 0.05 and 7 = 2), and fast parametrisations (/3 = 0.02 and 7 = 3). 



of the magnetic field B with the baryons density p from the center to the surface of a star is parametrized [14, 20] by 
the foUowing form 



Po 



1 — exp < — /? 



% 7 

Po 



(36) 



where po is the saturation density, B^^^^ is the magnetic field at the surface taken equal to lO^^G in accordance with 
the values inferred from observations and Bq represents the magnetic field for large densities. The parameters /3 and 
7 are chosen in such way that the field decreases fast or slow with the density from the centre to the surface. In this 
work, we will use two sets of values: a slowly varying field with /3 = 0.05 and j — 2, and a fast varying one defined by 
(3 = 0.02 and 7 = 3. We give the magnetic field in units of the critical field B^ = 4.414 x 10^^ G, so that B = B* B^. 
We further take Bq as a free parameter to check the effect of different fields. 

In Fig.l we plot the variation of the magnetic field, for the two parametrisations of the magnetic field as function 
of the baryonic density for several values of the parameter Bq in Eq. (36). For the slowly varying field the magnitude 
Bq is reached for densities p > lOpo, where po is the saturation density, while for the fast changing field this value 
is reached for p > 5po- On the other hand, for densities p < 3po the slowly varying magnetic field is stronger. This 
fact will influence the stiffness of the equation of state (EOS): the EOS is stiffer a) at low densities if a slowly varying 
field is considered; b) at high densities if a fast varying field is chosen. We present the numerical results for the 
slowly and fast parametrisations of the density-dependent magnetic field and for two Bag constants (ISOMeV)"' and 
(165MeV)^. In all figures where both cases are considered we show on the left panel the results for Bag^^^ —180 MeV 
and on the right panel the results for and Bag^/^=165 MeV. Furthermore, in all the figures, we will only show results 
obtained without the inclusion of the baryonic AMM's, because for the intensity of the magnetic fields considered the 
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FIG. 2: (Color online). EOS for stellar matter for GMl and for several values of the magnetic field identified by the Bq 
parameter of Eq. (36). In the left panel (a), (c), and (e) for Bag^'''* = 180 MeV and in the right panel (b), (d), and (f) for 
Bag^/*=165 MeV. In each figure we include the 5 = EOS, for reference, and the EOS obtained with the slow and fast 
parametrisation of Eq. (36). The hadron, mixed and quark phases are identified, respectevely, by a H, a M and thick lines, 
and a Q. 

contribution to the EOS is negligible. The effect of the AMM is non negligible for B* > 10^ [6] but, as it will be 
shown later in this section at the centre of the stars described in the present work these values will never be reached. 
In Fig. 2, we plot the EOS of neutron stellar matter, for the two parametrisations of the magnetic field, defined by 
Eqs. (28). For comparison we include in all figures the EOS for B = 0. The hadron phase is identified with a H 
(left part of the curves), the mixed phase is plotted with thick lines and identified with a M and the quark phase is 
identified with a Q (right part of the curves). We point out that these curves include both the contribution of stellar 
matter and the contribution of the magnetic field in the energy density and the pressure. In Fig. 3 wc plot only the 
pressure term coming from the stellar matter contribution, Eq. (27), as a function of the baryonic density for the 
largest Bq = 2 x 10^ in order to clearly show the effect of the magnetic field on the stellar matter contribution to the 
total EOS. We also plot, seperatly, the contribution of the magnetic field (thin line), which allows us to identify at 
which density the magnetic field contribution becomes larger than the matter contribution. The vertical lines identify 
the central density of the maxinmum mass star for both parameters bag considered. Figs. 2 and 3 are completed 
with Tables III and IV, where for the two parametrisations of the magnetic fields and, respectively, for Bag^/^=180 
MeV and Bag^/'*=165 MeV, the mixed phase density boundaries and the corresponding total energy densities defined 
in Eq. (28) are given for several values of the magnetic field. The onset of the mixed phase occurs at density pi, 
and the pure quark phase begins at density p2- In Tables V and VI we give, for comparison, the same quantities 
obtained using the Maxwell construction. Independently of the magnetic field, the onset of the mixed phase and 
quark pure phase occurs at lower densities for the smaller value of bag pressure [21], because the EOS is much stiffer 
for Bagi/'' = 180 MeV than Bag^^^ = 165 MeV. 

We now discuss the effect of the presence of a strong magnetic field with magnitude Bq = 5 x 10'*, 10^, and 2 x lO'"^ 
(~ 10^^ — 10^^ G) on the equation of state (EOS). The main conclusions are: a) the presence of the field does not 
affect the baryonic density at the onset of the mixed phase because the intensity of the field is too low (B* < 2 x ID"*); 
b) the baryonic density at the onset of the pure quark matter decreases with the increase of the magnetic field, as seen 
in Fig. 3 and discussed in Ref. [14]. This is due to the fact that the EOS of quark matter under a strong magnetic 
field becomes softer. The softening of the quark EOS with the increase of the magnetic field has been obtained with 
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p(fm"^) P(fm"^) 

FIG. 3: (Color online). EOS for stellar matter for GMl and for Bq ^ 2 x 10^ In the left panel (a) for Bagi/*=180 MeV and in 
the right panel (b) for Bag^/'*=165 MeV. In each fi gure we include the 5 = EOS, for reference, and the EOS obtained with 
the slow and fast parametrisation of Eq. (36). The vertical lines identify the central baryonic density of the stable star with 
maximum mass. 



TABLE III: The phase density boundaries, the onset of the mixed phase ui — pi/ po and the onset of the pure quark phase 
U2 ~ P2/P0, and the corresponding total energy densities, for several values of magnetic field using the two parametrisations 
defined in Eq. (36) (slow and fast). The bag constant is Bag^/* = 180 MeV, and the nuclear matter saturation density 
po = 0.153 fm"^ for GMl model. 



So 



Ml = Pl/po 



£i(fm 



"2 = P2/P0 



£2(fn 



B = 
BS = 5x 10* 

Bo* = 10' 
Be, = 2 X lO'' 



slow 
fast 
slow 
fast 
slow 
fast 



1.605 
1.604 
1.606 
1.609 
1.605 
1.614 
1.603 



1.212 
1.220 
1.218 
1.252 
1.228 
1.366 
1.273 



5.618 
5.588 
5.583 
5.510 
5.477 
5.357 
5.206 



4.794 
5.145 
5.333 
6.175 
6.914 
10.198 
13.017 



other quark models like the Nambu-Jona Lasinio (NJL) model [12, 22]; c) due to the magnetic field contribution 
to the pressure and energy density, Eqs. (28) , the extension of the mixed phase in terms of the energy density 
increases, as can be seen from Tables III and IV and from Fig. 2, mainly, shifting the onset of the quark phase to 
higher energy densities; d) the contribution of the magnetic field makes the total EOS harder. However, the slow 
and fast parametrisations have different effects. At low densities the soft parametrisation correponds to a stronger 
magnetic field and therefore gives rise to a harder EOS in that range of densities and favours the onset of the mixed 
phase at larger energy densities. In Tables III and IV Si is larger for the slow paramerisation. On the other hand, 
for large densities the fast parametrisation corresponds to larger magnetic fields and a harder EOS. The onset of the 
quark phase occurs at larger energy densities, i.e. for a given Bq value, £2 is larger for the fast parametrisation; e) 
the transition to the mixed phase is less affected by the magnetic field than the transtion to a pure quark matter. 
This is due to the magnetic field parametrisation, (36), which gives a weaker magnetic field at the baryonic densities 
of the mixed phase onset. On the other hand, the transition to a pure quark phase is strongly affected because for 
the densities at which it occurs the magnetic field is stronger, more than 50% (75%) its maximum value for the slow 
(fast) varying magnetic field; f) for all values of the magnetic field considered, the extension of the mixed-phase is 
larger for the larger bag pressure due to the larger stiffness of the quark phase EOS. From Tables V and VI, for the 
Maxwell construction similar conclusions are drawn, however with a larger baryonic density at the onset of the mixed 
phase, ^ Po larger for Bag^/^=180 MeV and ~ 0.5po larger for Bag"'^/'*=165 MeV, and a smaller one at the on set of 
the pure quark phase, ~ po smaller for Bag^/''=180 MeV and ~ O.Spo smaller for Bag^/'*=165 MeV. 
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In Rcf. [14] the magnetic field contribution was also taken into account in the EOS used to integrate the TOV 
equations, however this term enters with the wrong sign in the total pressure, it gives a negative contribution, and 
therefore, the EOS becomes very soft for strong magnetic fields and the formation of very massive stars is not allowed. 



TABLE IV: Same as Table III, but for Bag^/*=165 MeV. 

Bg m = Pi/po £i(fm~'^) U2 = p-z/po £2(fm~^) 

B = 0.832 0.613 3.567 2.748 

Bo = 5 X 10* slow 0.831 0.613 3.546 2.869 

fast 0.831 0.612 3.538 2.940 

So = 10^ slow 0.830 0.614 3.519 3.232 

fast 0.831 0.612 3.478 3.467 

= 2 X IQp slow 0.830 0.623 3.402 4.496 

fast 0.832 0.614 3.317 5.162 



TABLE V: The phase density boundaries obtained with Maxwell construction, the onset of the mixed phase wi = pi/po and 
the onset of the pure quark phase U2 = P2/P0, and the corresponding total energy densities, for several values of magnetic field 
using the two parametrisations defined in Eq. (36) (slow and fast). The bag constant is Bag^'''*=180 MeV, and the nuclear 
matter saturation density po = 0.153 fm~"^ for GMl model. 



Ml 



■ Pl/PO 



ei(fm-4) 



U2 



P2/P0 



e2(fm 



5 = 
S5 = 5 X 10"* 

So* = 10^ 

Bo = 2 X 10^ 



slow 
fast 
slow 
fast 
slow 
fast 



2.819 
2.799 
2.787 
2.706 
2.715 
2.477 
2.475 



2.250 
2.297 
2.297 
2.379 
2.422 
2.628 
2.612 



4.606 
4.580 
4.580 
4.510 
4.471 
4.228 
4.314 



3.855 
4.089 
4.275 
4.760 
5.420 
6.918 
9.813 



TABLE VI: Same as Table V, but for Bagi/''=165 MeV. 

Bq Ml = pi/po £i(fm~"^) U2 = P2/P0 £2(fm~'^) 

B = 1.283 0.958 2.922 2.235 

B5 = 5 X 10* slow 1.276 0.956 2.933 2.318 

fast 1.276 0.953 2.933 2.339 

Bo = 10^ slow 1.260 0.954 2.918 2.524 

fast 1.255 0.939 2.902 2.583 

Bo = 2 X 10^ slow 1.206 0.946 2.878 3.321 

fast 1.196 0.901 2.847 3.507 



As already discussed in Ref. [21], the presence of strangeness in the core of a neutron stars will affect some of their 
properties. In Fig. 4 we show the strangeness fraction defined by : 

t. = X^T + (1 - X)tr (37) 

with, 

= = ^^M^ (38) 

where pa is the strange quark density, and q\ is the strange charge of baryon 6, listed in Table I. 

Only the strongest field considered affects considerably the strangeness fraction. At the onset of the pure quark 
phase, the strange quark fraction has reached 30% of the baryonic matter. Although the fraction of strangeness varies 
in the mixed-phase, it has almost the same value in the quark phase. For Bq = 2 x 10^, with both bag values and with 
the two parametrisations of the magnetic field used the effect of the Landau quantisation appears, and the strangeness 
fraction obtained is slightly oscillating around the B = Q result. The strangeness content increases slightly faster at 
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larger densities with the fast parametrisation. The presence of this strong magnetic field shifts to lower densities the 
onset of the 30% fraction of strange matter, corresponding to the onset of the quark phase. The thin vertical lines 
included in the Fig. 4 identifies the central baryonic density of the maximum mass star. For Bag^/*=180 MeV the 
strangeness fraction at given density is never larger than 15%, while for Bag-'^/'*=165 MeV it can be as large as 30%. 



TABLE VII: Properties of the stable hybrid star with maximum mass, for Bag^/* = 180 MeV and for several values of magnetic 
field using [see Eq. (36)] with the slow and fast parametrisations. Mmax, M^ax, R, Eo, p'^, and B* are, respectively, the 
gravitational and baryonic masses, the star radius, the central energy density, the central baryonic density, and the value of 
the magnetic field at the centre. 



Bo* 




Mmax[MQ\ 




R [km] 


£;o[fm-''] 


= p^/po 


Bl 


B = 




1.42 


1.57 


10.22 


8.05 


8.680 




55 = 5 X 10'' 


slow 


1.74 


1.93 


11.56 


5.75 


6.150 


4.248 xlO* 




fast 


1.85 


2.06 


11.49 


5.85 


6.092 


4.948 xlO"* 


Bo* = 10' 


slow 


2.23 


2.48 


12.35 


5.23 


4.843 


6.908 xlO* 




fast 


2.33 


2.55 


11.88 


5.62 


4.595 


8.565 xlO'' 


B5 = 2 X lO'' 


slow 


2.86 


3.12 


14.39 


3.73 


3.033 


7.375 xlO" 




fast 


2.76 


2.88 


13.15 


4.41 


3.092 


8.926 xlQ'' 



TABLE VIIL Same as Table VII, but for Bagi/''=165 MeV. 



Bo* 




Mrr,ax[MQ\ 


MlaxWe] 


R [km] 


£o[fm-^] 


u" = p7po 


Bt 


B = 




1.47 


1.68 


8.86 


10.02 


10.673 




B5 = 5 X lO-* 


slow 


1.72 


1.95 


9.80 


8.24 


8.588 


4.877 xlO"* 




fast 


1.81 


2.04 


10.07 


7.55 


7.928 


5.000 xlO* 


BS = 10' 


slow 


2.14 


2.36 


11.07 


6.52 


5.902 


8.251 xlO"* 




fast 


2.26 


2.45 


11.01 


6.36 


5.203 


9.404 xlO'' 


Bo = 2 X lO'' 


slow 


2.73 


2.93 


13.55 


4.23 


3.294 


8.377 xlO"* 




fast 


2.68 


2.78 


12.50 


4.80 


3.222 


9.759 xlO* 



We next study the effect of the magnetic field on the properties of the stars described by the EOS discussed above. 
We will consider that the stars have spherical symmetry. The maximum gravitational and baryonic masses of the 




FIG. 4: (Color online). Strangeness fraction, defined in Eq. (37), as function of the baryonic density, for Bq = 2 x lO' in 
the interior of the star, (a) for Bag'^^'' = 180 MeV and (b) for Bag'^/'' = 165 MeV. The vertical lines identify the central baryonic 
density of the stable star with maximum mass. 
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stable stars, radius, central energy density, and central baryonic density, and the corresponding magnetic field at the 
centre of the star, obtained by solving the Tolman-Oppenheimer-VolkofF equation, are given for Bagi/4=180 MeV and 
Bag^''^=165 MeV in Tables VII and VIII, respectively. We note that the central energy and baryonic densities are 
decreasing when the magnetic field increases, while the maximum masses and the corresponding radius increase with 
the magnetic field, contrary to the results of [14] where the magnetic field contribution has entered with the wrong 
sign in the total pressure expression, giving rise to very soft EOS for large fields. This behavior is due to the stiffening 
of the EOS with the magnetic field. In the cases under study the magnetic field at the center of the maximum mass 
configuration is never larger than 4 x 10^* G a value close to the limit obtained from the virial theorem. 



— B=0 




R [km] R [km] 

FIG. 5: (Color online). Mass-radius relation of the hybrid stars described in the present work for several values of magnetic 
field using the slow and fast varying parametrisations of B. (a) and (c) for Bag'^/*=180 MeV, and (b) and (d) for Bag^^'* = 165 
MeV. The thick lines identify the stars with a mixed phase at the center. Stars with smaller masses are hadronic stars with no 
quark matter, and stars with larger masses have a pure quark core. For the Bq = 2 x lO"" there are no stars with a quark core. 

The mass/radius curves for the families of stars corresponding to the maximum mass configurations given in 
Tables VII and VIII are plotted in Fig. 5. The thick lines correspond to stars for which the central density lies 
within the mixed phase. For the largest bag pressure there is a quark core only for the no field configuration and the 
smallest field considered. For the smallest bag pressure only the strongest field considered gives rise to a no quark core 
configuration. It is interesting to notice the strong effect of the magnetic field parametrisation on the corresponding 
family star properties. For the stars that have no quark phase in their interior the central density is always smaller 
than p ~ 2.5po for which the slow and the fast varying parametrisations give the same magnetic field magnitude, 
see Fig. 1. Below this density the slowly varying parametrisation gives rise to larger magnetic fields and therefore 
to a harder EOS. The corresponding stars have larger radius and larger masses due to the larger incompressibility 
modulus. This explains why the maximum radius of stars with more than 0.5 Mq is larger for the slowly varying 
field. This difference can be as large as 1 km for the largest field considered. 

However, the largest mass configuration is obtained for central magnetic fields larger for the fast varying field and 
therefore, the maximum mass is larger in this case, except for the largest central value of the magnetic field considered. 
For this value, the central density for the maximum mass configuration is 3 — 3.3 po- From Fig. 1 it is clear that 
these densities are just above the critical density values for which both parametizations give the same magnetic field 
and, therefore, the star obtained with the slowly varying parametrisation has, for most of the densities a larger field. 
The size of the star is largely influenced by the lower density layers and therefore most of the maximum mass stars 
have a larger radius for the slowly varying field, which give rise to larger magnetic fields in the low density layers. 

For a central magnetic field ~ 3 x 10^^ G we get maximum mass configuration with a mass 2.2-2.3 Mq. Slightly 
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larger central fields predict even more massive stars with M > 2.7 Mq. These values would be able to describe highly 
massive compact stars, such as the one associated to the millisecond pulsars PSR B1516 + 02B [23], and the one 
in PSR J1748-2021B [24] in case they are confirmed. Otherwise an upper limit on the possible magnitude of the 
magnetic field at the center of a compact star may be obtained. 




FIG. 6: (Color online). Particle fractions Yi — pi/p as function of p/po, for i=baryons, leptons, and quarks, for two values of 
the magnetic field (Bo = 5 x 10'' (a)-(d), and Bq — 2 x 10^ (e)-(h)), with both parametrisations (slow and fast). We have 
considered two bag pressures: Bag^''* = 180 MeV for (a), (c), (e), and (g) figures, and Bag''^'=165 MeV for (b), (d), (f), and 
(h) figures. The vertical lines represent the density p'^ at the center of the star with majcimum mass. 



The properties of the stars are strongly influenced by the magnitude of the magnetic field in their interior and the 
characteristics of the hadronic EOS become less important the larger the field is. It is also clear that a quark phase 
is not favoured if a strong magnetic field exists in the interior of the star. This conclusion is based on the Gibbs 
construction of the mixed phase. Within a Maxwell construction, the onset of pure quark matter occurs at smaller 
densities and, therefore, a small quark core could still exist for the smaller Bag constant. 

We have also looked at the baryonic and leptonic composition of the stars. In Fig. 6 we plot the lepton, baryon and 
the quark abundances. We only show the particle fractions for the smallest and the largest magnetic fields considered 
in the present work. The vertical lines in the figures represent the central density of the configuration with maximum 
mass when it lies within the range of densities shown. Large magnetic fields prevent the appearance of a quark phase, 
enhance the leptonic fraction and reduce slightly the mixed phase baryonic density extention. For strong magnetic 
fields, stars are mainly hadronic even when we consider a soft quark EOS. We also conclude that the only effect of the 
slow/fast varying magnetic field is the appearance at low/high densities of a larger leptonic fraction. The contribution 
of hyperons is minimal or even zero. 



VI. CONCLUSIONS 



We have studied the influence of a static magnetic field on the deconfinement phase transition in the interior of a 
compact star. For the hadronic phase we took the GMl parametrisation of the NLWM [16] and for the quark phase 
the MIT bag model with two different bag pressures, one corresponding to a soft and the other to a hard quark EOS. 
The mixed phase was built imposing the Gibbs conditions [13]. It was considered that the magnitude of the magnetic 
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field increases with density from a surface field of the order of 10^^ G to a central field ^ 5 x lO^^G, and both a fast 
and a slow changing field were used. 

It was shown that due to the weakness of the magnetic field at low densities, the onset of the mixed phase is not 
affected by the magnetic field. The baryonic density at the onset of the pure quark phase is only slightly reduced 
for the magnetic fields considered, but, due to the contribution of the magnetic field to the total pressure and total 
energy density, the EOS of stellar matter in the presence of a strong magnetic field becomes much harder and the 
mixed phases extends to larger energy density ranges. As a consequence, the family of stars corresponding to these 
EOS have configurations with larger masses and radius. It was shown that the density dependence of the magnetic 
field will influence the star properties, with a slowly increasing field giving rise to larger radius. 

Due to the magnetic field pressure, the quark phase is not favoured, and stars with strong magnetic field are mainly 
hadronic with very small hyperon contributions. Stars with very high masses and radius are predicted and maximum 
masses of observed compact stars may set an upper limit of the largest possible magnetic field at the center of the 
star, - 2 X 10^* G for 2 Mq stars. 

In the present work we have not taken into account the possibility of color superconductivity (CS) which was shown 
to be important to describe the groundstate of quark matter at high chemical potential and low temperature (for 
a review see [25]). In particular, the effect of CS, in a color-flavor locked (CFL) phase, on the structure of hybrid 
stars was discussed in Ref. [26]. Moreover, it has already been studied the effect of strong magnetic fields on the 
CS properties which can be drastic, [27]. However, we point out that although a strong magnetic field may have 
important effects on the stellar quark matter EOS the contribution of the magnetic field energy to the total energy of 
the star will impose severe restrictions on the maximum field a star may support and on the baryonic density reached 
at the centre of the star, [15]. In the present work, we have shown that: for very strong magnetic fields the central 
core of the star will probably never be in a quark phase, but at most in a mixed phase. 

Very strong magnetic fields can only occur in very young compact stars before the magnetic field has decayed [28]. 
We may therefore conclude that, not only neutrino trapping in the protoneutron star phase [29], but also a very strong 
magnetic field in the interior of a compact star strongly precludes the formation of a quark phase. Quark matter 
thermal nucleation in hot and dense hadronic matter has been proposed by several authors [30]. In these studies, it 
was found that the prompt formation of a critical size drop of quark matter via thermal activation is possible above 
a temperature of about 2 — 3 MeV, and as a consequence, it was inferred that pure hadronic stars are converted to 
hybrid stars or quark stars within the first seconds after their birth. The presence of a strong magnetic field would 
make this scenario less probable. 

Acknowledgments 

This work was partially supported by FEDER and Projects PTDC/FP/64707/2006 and CERN/FP/83505/2008, 
and by COMPSTAR, an ESF Research Networking Programme. 



[1] R. C. Duncan and C. Thompson, Astronphys. J. 392, L9 (1992). 

[2] V. V. Usov, Nature 357, 472 (1992). 

[3] B. Paczynski, Acta Astron. 42, 145 (1992). 

[4] S. L. Shapiro and S. A. Tcukolsky, Black holes, white dwarfs and neutron stars, Wiley-interscience New York, 1983. 
[5] S. Chakrabarty, Phys. Rev. D 54, 1306 (1996). 

[6] A. Broderick, M. Prakash, and J. M. Lattimer, Astrophys. J. 537, 351 (2000). 

[7] J. Boguta and A. R. Bodmer, Nucl. Phys. A292, 413 (1977); B. D. Serot and J.D. Walecka, Adv. Nucl. Phys. 16, 1 (1986). 
[8] A. Rabhi, C Providencia, and J. da Providencia, J. Phys. G: Nucl. Part. Phys. 35,125201 (2008). 
[9] F.X. Woi, G.J. Mao, CM. Ko, L.S. Kisslingcr, H. Stocker and W. Greiner, J. Phys. G 32, 47 (2006). 

[10] H. Pais, Master thesis, University of Coimbra, 2008. 

[11] R. Gonzalez Felipe, A. Perez Martinez, H. Perez Rojas and M. Orsaria, Phys. Rev. C 77, 015807 (2008). 

[12] D. P. Menezes, M. Benghi Pinto, S. S. Avancini, A. Perez Martinez, and C. Providencia, Phys. Rev. C 79, 035807 (2009). 

[13] N. K. Glendenning, Compact stars, Springer- Verlag New York , 2000. 

[14] D. Bandyopadhyay, S. Chakrabarty, and S. Pal, Phys. Rev. Lett. 79, 2176 (1997). 

[15] C. Y. CardaU, M. Prakash, and J. M. Lattimer, Astrophys. J. 554, 322 (2001); A. Broderick, M. Prakash and J.M. Lattimer, 

Phys. Lett. B 531, 167 (2002). 
[16] N. K. Glendenning and S. A. Moszkowski, Phys. Rev Lett. 67, 2414 (1991). 

[17] D. N. Voskresensky, M. Yasuhira, and T. Tatsumi, Phys. Lett. B 541, 93 (2002); D. N. Voskresensky, M. Yasuhira, and T. 

Tatsumi, Nucl. Phys. A723, 291 (2003). 
[18] T. Mamyama, S. Chiba, H-J Schulze, and T. Tatsumi, Phys. Rev. D 76, 123015 (2007). 



15 



[19] M. Bocquet, S. Bonazzola, E. Gourgoulhon, and J. Novak, Astron. Astrophys. 301, 757 (1995). 

[20] Guang-Jun Mao, Akira Iwamoto, and Zhu-Xia Li, Chin. J. Astron. Astrophys. Vol 3, No. 4, 359-374 (2003). 

[21] D. P. Menezes and C. Providencia, Phys. Rev. C 68, 035804 (2003); I. Bombaci, I. Parenti, I. Vidana, Astrophys. Jour. 

614, 314 (2004). 
[22] D. Ebert, K.G. Klimenko, Nucl. Phys. A 728, 203 (2003). 

[23] Paulo C. C. Preire, Alex Wolszczan, Maureen van don Berg, and Jason W. T. Hessels, Astrophys. J. 679, 1433 (2008). 
[24] Paulo C. C. Preire, Scott M. Ransom, Steve Bgin, Ingrid H. Stairs, Jason W. T. Hessels, Lucille Prey, Fernando Camilo, 

Astrophys. J. 675, 670 (2008). 

[25] D. Bailin and A. Love, Phys. Rep. 107, 325 (1984); K. Rajagopal and F. Wilczok, arXiv:hep-ph/0011333v2; M. Alford, 
Ann. Rev . Nucl. Part. Sci. 51, 131 (2001); Riv . Nuovo Cim. 25N3, 1 (2002); T. Schaefer, arXiv:hep-ph/0304281v2; D.H. 
Rischke, Prog. Part. Nucl. Phys. 52, 197 (2004); Hai-cang Ren, arXiv:hep-ph/0404074v3. 

[26] M. Alford and S. Reddy, Phys. Rev. D 67, 074024 (2003); P. K. Panda, D. P. Menezes, and C. Providencia, Phys. Rev. C 
69, 025207 (2004); P. K. Panda, D. P. Menezes, and C. Providencia, Phys. Rev. C 69, 058801 (2004). 

[27] K. lida and G. Baym, Phys. Rev. D 66, 014015 (2002); E. J. Ferrer, V. de la Incera and C. Manuel, Phys. Rev. Lett. 95, 
152002 (2005); E. J. Ferrer, V. dc la Inccra and C. Manuel, Nucl.Phys. B 747, 88 (2006); E. J. Ferrer and V. do la Inccra, 
Phys. Rev. Lett. 97, 122301 (2006); E. J. Ferrer and V. de la Incera, Phys. Rev. D 76, 045011 (2007); E. J. Ferrer and 
V. de la Incera, Phys. Rev. D 76, 114012 (2007) ; J. L. Noronha and I. A. Shovkovy, Phys. Rev. D 76, 105030 (2007); K. 
Fukushima and H. J. Warringa, Phys. Rev .Lett. 100, 032007 (2008); D. M. Sedrakian, K. M. Shahabasyan, D. Blaschke 
and M. K. Shahabasyan, Astrophysics, Vol. 51, No. 4, (2008). 

[28] http:/ /www. physics. mcgill.ca/~pulsar/magnetar/main. html 

[29] I. Bombaci, I. Parenti and I. Vidana, Astrophys. J. 614, 314 (2004); Ignazio Bombaci, PrafuUa K. Panda, Constanca 

Providencia, and Isaac Vidana, Phys. Rev. D 77, 083002 (2008). 
[30] J. E. Horvath, O. G. Benvenuto and H. Vucetich, Phys. Rev. D 45, 3865 (1992); J. E. Horvath, Phys. Rev. D 49, 5590 

(1994); M. L. Olesen and J. Madsen, Phys. Rev. D 49, 2698 (1994). 



